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Background Nuclear fusion has been shown to be a useful probe to study the different nuclear shapes. However, the pos¬ 
sibility of testing octupole deformation of a nucleus with this tool has not been fully explored yet. The presence of a 
stactic octupole deformation in nuclei will enhanced a possible permanent electric dipole moment, leading to a possible 
demonstration of parity violation. 

Purpose To check whether static octupole deformation or octupole vibration in fusion give different results so that both 
situations could be experimentally disentangled. 

Method Fusion cross sections are computed in the Coupled-Channels formalism making use of the Ingoing-Wave Boundary 
Conditions (IWBC) for the systems 16 0+ 144 Ba and 16 0+ 224 Ra. 

Results Barrier distributions of the two considered schemes show different patterns. For the 224 Ra case, larger differences are 
found in correspondence with its larger octupole deformation. 

Conclusions The measurement of barrier distributions can be a complementary probe to clarify the presence of octupole 
deformation. 


PACS numbers: 25.60.Pj, 24.10.Eq 


I. INTRODUCTION 

The search of octupole deformations in nuclei is liv¬ 
ing a revival thanks to its impact in a possible perma¬ 
nent atomic electric dipole moment (EDM) 0,1. A 
non-zero EDM will indicate a time-reversal (or equiva¬ 
lently charge-parity) violation. The magnitude or the 
experimental maximum limit to it can constraint the dif¬ 
ferent suggested extensions to the standard model Q. 
The presence of a static octupole deformation in an odd 
nucleus will generate enhanced nuclear Schiff moments, 
which contributes to the atomic EDM so that it can be 
improved by several orders of magnitude 0]. 

Therefore, the experimental focus is set on the search 
of permanent octupole deformations in some of the 
regions where different theoretical approaches predict 
strong octupole correlations 0, [|] with the help of con¬ 
tinuous development of the radioactive beam facilities. 
In the present work we will focus on 144 Ba and 224 Ra 
as representative of two different regions where possible 
static deformations have been predicted. 

In 0,0 , Coulomb excitation has been used to mea¬ 
sure the different electric transition probabilities. This 
tool provides quadrupole and octupole transition proba¬ 
bilities with good accuracy. Smaller dipole electric transi¬ 
tion probabilities will carry larger error bars, even though 
this problem is experimentally affordable and they were 
able to provide some measurements. It should be kept in 
mind that large octupole transitions can be found also for 
dynamic octupole vibrations. In addition, the coupling 
between quadrupole deformation and octupole vibrations 
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can lead to enhanced dipole moments 0. 

Therefore, we would like to propose here a comple¬ 
mentary experimental probe for static octupole deforma¬ 
tion. A traditional experimental tool for the study of 
nuclear structure is provided by subbarrier fusion. It is 
well known that fusion at energies around the Coulomb 
barrier is driven by the dynamical couplings to the inter¬ 
nal degrees of freedom of the fusing counterparts @, 0. 
The absence or presence of octupole and dipole moments 
in one of the fusing partners will have a certain impact in 
the final subbarrier fusion cross section. This fact could 
suggest, if the minimum accuracy is reached, the possibil¬ 
ity of distinguishing between static octupole deformation 
and the corresponding dynamical vibration. 

The work is structured as follows. In Sec. El we re¬ 
call the reaction formalism for studying nuclear fusion 
including deformations and/or vibrations. We apply 
this framework in Sec. IIIII to the reactions 16 0+ 144 Ba 
and 16 0 + 224 Ra considering quadrupole deformations for 
144 Ba and 224 Ra looking for the differences between 
adding an octupole vibration or deformation to the pre¬ 
vious quadrupole deformation. Finally, in Sec. ITVl the 
main results of this work are summarized. 


II. REACTION FRAMEWORK 

Fusion probabilities are calculated by solving the corre¬ 
sponding coupled-channel equations under ingoing-wave 
boundary conditions (IWBC). The coupled-channel for¬ 
malism for direct reaction processes given by Austern m 
expands the total wave function in terms of the wave- 
function for the internal state of the projectile and 
the radial wave functions XB that accounts for the rela¬ 
tive motion between projectile and target. This leads to 
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TABLE I: Theoretical deformation parameters for 144 Ba and 
224 Ra according to [l3| for the different multipolarities. 



144 Ba 

224 Ra 

fa 

0.149 

0.138 

P3 

0.068 

0.099 


a set of coupled equations for the radial wave functions: 


d 2 Xp 

dR 2 


L(L + 1) 2/i 

' U H 2 


R 2 


+ 3W - vf’mxt = 
= % E V^’(R)x. 

CL^fi 


(1) 


In these expression V is the interaction potential, p is the 
reduced mass, and Ep is the relative energy. Each chan¬ 
nel /3 corresponds to a set of quantum numbers {/, L, J}, 
where I is the angular momentum of the internal state, 
L is the angular momentum in the relative coordinate, 
and both are coupled to a total angular momentum J. 

As a simplification we can use the same central poten¬ 
tial for the (R) although this potential is modified 
by the terms Vpp known as reorientation terms. The non 
central terms Vp a are those on charge of the coupling be¬ 
tween two different channels a and /3. For all the cases, 
we will consider, for the central potential, a Woods-Saxon 
with the parametrization from Akyiiz-Winther EMI 
plus the corresponding Coulomb repulsion. 

In our case we will consider one of the two ions to be 
a quadrupole deformed rotor.We will also look at cases 
where the deformed nuclei has an octupole vibration or a 
static octupole deformation together with the quadrupole 
one. 

In the case of permanent deformation, we can describe 
the radius of the deformed nucleus as a function of the 
angle 6', defined with respect to an intrinsic body-fixed 
frame, R{0') = i?o[l + PxY£ 0 (9')\, where R 0 is an 

average radius of the deformed nucleus and /3\ a dimen¬ 
sionless deformation parameter. 

The nuclei of interest in the present work are candi¬ 
dates to have a permanent octupole deformation together 
with the quadrupole one. In Table Q] we show the defor¬ 
mation parameters used in this work following the theo¬ 
retical predictions from (Td • 

If one assumes that the potential is still a function of 
the distance between projectile and target, the potential 
can be expanded in multipoles as: 

v(R,n) = Y / v x(R)' D U a ’P^) Y x l .(R), ( 2 ) 

A/i, 

where T> is the so called rotation matrix (or D- 
matrix) 0- Finally, evaluating the matrix elements of 
this potential between the states of the rotor, it is pos¬ 
sible to obtain the coupling potentials. See @,10,0 
for more details. 


The spherical harmonic will connect two states x a (ft) 
and X/?(i?) depending on the order A. The iso-centrifugal 
approximation ilGlll8l | is used to reduce the size of the 
calculation. 

On the other hand, the coupling may also arise from 
a vibration of one of the nucleus. In this case this vi¬ 
bration is characterized as a variation on the surface 
as i?(£) = i? 0 [l + where a X p are to 

be understood as dynamical variables, given in terms 
of phonon creation (6^ ) and annihilation (fix^) opera¬ 
tors 0. The nuclear coupling between the ground state 
and the first one phonon state of multipolarity A reads: 

Vcoup = (3) 

However, we will consider an octupole vibration on top of 
a quadrupole deformed nucleus instead of a pure spher¬ 
ical nucleus. In this case, the derivative of the potential 
has a certain dependence of the orientation. We can ex¬ 
press such dependence in the form: 

A, dV 

Vcoup(R,0 = — (ft — fto [1 + /32^2 o]) ^30, (4) 

which we later expand in spherical harmonics as done 
for the full rotor case, i.e. following Eq. ©• Here, 
we have selected p = 0 because the coupling of the oc¬ 
tupole vibration with the quadrupole deformation splits 
the strength into different bands, being K = 0” the low¬ 
est one 0, 0]. In other words, if the nuclei is signifi¬ 
cantly prolate, the vibration is stronger along the sym¬ 
metry axis where the radius is larger. One should remind 
that this potential is only coupling zero octupole phonon 
states with one octupole phonon states. Together with 
these vibrational couplings, the traditional quadrupole 
deformed potential will only act between states with the 
same number of phonons. 

In addition to all these nuclear couplings considered 
here, high order Coulomb couplings should be consid¬ 
ered. However, it has been shown that Coulomb cou¬ 
plings has a minor role compare to the nuclear ones [22j 
and, therefore, they will not be included here. 

Finally, solving the coupled equations a probability of 
transmission Tl for each angular momentum L is ob¬ 
tained. The fusion cross-section is given by 

(5) 

L=0 " L=0 


The probability of transmission for each partial wave can 
also be calculated simply by a shift of energy, 


Tr 5* Tn 


E- 


L(L + l)h 2 
2pR^ 


( 6 ) 


where Rq is the position of the barrier for the s-wave [8}. 
Hereafter we will refer as transmission probability T to 
the probability of transmission for the s-wave To jlj, 0] . 













3 




E (MeV) 

c.m. v ' 


FIG. 1: (Color online) Barrier distribution for the fusion of 
16 O and 144 Ba including a positive quadrupole deformation 
for 144 Ba, solid line, and the same for a negative quadruople 
deformation, dashed line. 


The derivative of the transmission probability, usually 
called barrier distribution, is approximately proportional 
to the second derivative of the product of the cross section 
and the energy ( Ea ), thus being a direct link between nu¬ 
clear structure and this experimental observable mu¬ 


ni. APPLICATIONS 


A. Prolate vs. Oblate deformations 



(MeV) 


FIG. 2: (Color online) Cross sections (upper panel) for the 
fusion of 16 0 and 144 Ba considering 144 Ba as a spherical nu¬ 
clei (solid line) and including a positive quadrupole deforma¬ 
tion within the frozen approximation (dashed-dotted line) and 
within the coupled channels framework (dashed line). For the 
last two, barrier distributions are shown in the lower panel. 


Fusion reaction is one of the experimental tests that 
can discriminate between prolate and oblate deforma¬ 
tions. A change on the sign of the deformation length 
will change the sign of the coupling and reorientation 
terms. The sign of the off-diagonal coupling potential 
does not have an effect in the cross section. Instead, the 
reorientation term is on-diagonal, so being the main re¬ 
sponsible for the change in the cross section when one of 
the counter partners is either oblate or prolate. 

In Fig. |T] we show the barrier distribution for the 
fusion of the system 16 0+ 144 Ba considering a positive 
quadrupole deformation parameter, solid line, and a neg¬ 
ative one, dashed line. In this calculation and here- 
afer, we consider only the first three leves of a typical 
quadrupole deformed rotor, including all allowed cou¬ 
plings. The influence of higher energy and spin levels 
depends on the nucleus. For 144 Ba we have studied the 
variation of the barrier distribution with the number of 
levels included. From this analysis we concluded that in¬ 
cluding more levels does not alter dramatically the bar¬ 
rier distribution. Similar results are found for 224 Ra. 

Additional examples of fusion reactions with prolate 
and oblate nuclei can be found in @, H3 


B. Coupled-Channels vs. frozen approximation 

The effect of static deformation was previously stud- 
ied in terms of the frozen approximation in [261 ] for 
quadrupole and in (27] for octupole deformation. Within 
this formalism the energies of the different states are ne¬ 
glected. As a consequence, the total cross section can 
be defined as the average cross section between those ob¬ 
tained in a single channel calculation for the different bar¬ 
riers found for each possible orientation of the deformed 
nucleus. This fact simplifies the calculations but it has 
been shown to overestimate the cross sections at energies 
below the barrier HU- However, the first states of the 
nuclei investigated have small excitation energies, what 
should reduce the possible overestimation. 

In order to test the validity of the frozen approxima¬ 
tion for these nuclei, we have studied the fusion of 16 0 
and 144 Ba considering only the quadrupole deformation. 
In the Fig. [2] we show the result of the Coupled Channel 
calculation as in the previous subsection and the result of 
considering the frozen approximation for the same defor¬ 
mation parameter fa ■ We see in Fig. UKa) that, at ener- 
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FIG. 3: (Color online) Level schemes considered for a 
quadrupole octupole deformed nucleus (left panel) and for a 
quadrupole deformed nuleus with an octupole vibration (right 
panel). Each arrow can be related to transitions with more 
than one possible multipolarity A. 



FIG. 4: (Color online) Coupling potentials for the system 
16 0+ 144 Ba considering quadrupole octupole deformed 144 Ba 
(left panel) and a quadrupole deformed nuleus with an oc¬ 
tupole vibration (right panel). 


gies below the barrier, the frozen approximation slightly 
overestimates the cross section. Barrier distributions are 
shown in Fig. EJb). Main qualitative features are kept so 
that the frozen approximation is consistent with the cou¬ 
pled channel calculations for this case. Therefore, this 
aproximation can still be of great interest but should be 
managed with certain caution. 


C. Difference between octupole deformation and 
vibration 

A more interesting and up-to-date case is to check if 
barrier distributions are sensitive enough to the difference 
between a nucleus with both quadrupole and octupole 
static deformations and a quadrupole deformed nucleus 
with an octupole vibration. The levels and the coupling 
potentials between the two cases will be different. In 
Fig. [3] we show the two different coupling schemes con¬ 
sidered. Each arrow represents the presence of at least 
one coupling potential with a certain multipolarity. 
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FIG. 5: (Color online) Barrier distribution for the fusion of 
16 0 on 144 Ba (upper panel) and 224 Ra (lower panel) including 
a quadrupole and octupole deformation for the target (solid 
line) or a quadrupole deformation plus an ocutpole vibration 
(dashed line). 


The scheme (a) in Fig. [3] consists of a single octupole- 
quadrupole deformed band. Blue solid arrows in this 
scheme are there to remind that we can couple all this 
levels with a deformed potential considering the nucleus 
both octupole and quadrupole deformed. The scheme 
(b) in Fig. [3] consists of a ground state quadrupole de¬ 
formed rotor band and a one octupole phonon quadrupole 
deformed band. Different one phonon bands may re¬ 
sult from the coupling of one octupole phonon and a 
quadrupole deformation. For big enough values of 02 
and & the lowest energy band is the K = 0~ band with 
associated levels I = 1 — , 3 — , 5 ... [2Cj |2l| so that this 
case coincides with the octupole quadrupole deformed 
level scheme [3, EH] • 

Only red dashed arrows connect these two bands since 
we need a one phonon creation to go from the ground 
state to any of the levels in this second band. 

The strength and shape of each coupling potential for 
the system 16 0+ 144 Ba for each multipolarity are shown 
in Fig. [4] again for the two different schemes considered. 
Differences between both schemes are small but signifi¬ 
cant. For a meaningful comparison the same value has 
been used for the static 03 in case (a) and the dynami¬ 
cal 03 in (b). The presence of an octupole deformation 
affects to the quadrupole strength and viceversa. Fig. [4] 
only includes up to A = 6. However, larger multipolari¬ 
ties have been included. 

Even though the energies expected for the negative 
parity states differ in each scheme, we will use in all cal¬ 
culations the experimental energies [30j collected in Ta¬ 
ble HD Doing so, we put ourselves in the less favourable 
situation since the difference in energy will translate into 
larger differences in the barrier distribution. 

With the procedure discussed in section ID we calcu- 
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TABLE II: Experimental excitation energies in MeV for 144 Ba 
and 224 Ra states according to [3Cj ]. 


r 

1" 

2+ 

3“ 

4+ 

5" 

6+ 

144 Ba 

0.759 

0.199 

0.838 

0.530 

1.038 

0.962 

224 Ra 

0.216 

0.084 

0.290 

0.250 

0.433 

0.479 


lated the barrier distributions for two systems 16 0+ 144 Ba 
and 16 0 + 224 Ra. Both 144 Ba and 224 Ra are candidates 
to have an octupole deformation with a considerable (3%. 
Excitations of the projectile are not considered. Results 
are shown in Fig. [31 Barrier distributions for both cases 
show some differences, but probably not enough to open 
the possibility of clearly distinguishing the two situations. 
It would depend on the strength of the octupole defor¬ 
mation and the experimental accuracy available on each 
case. However, it will be worthy to extend the present 
analysis to other regions with large octupole deforma¬ 
tions. As just an example, 232 Th could have an even 
larger =0.14 Hlj as investigated by K. Hagino and 
K. Sannohe f3l| . 

IV. CONCLUSIONS 

In this work we have studied the fusion reactions 
16 0+ 144 Ba and 16 0+ 224 Ra under different assumptions 
for the structure of 144 Ba and 224 Ra. We started intro¬ 
ducing only a quadrupole deformation and studying the 
sensitivity of the barrier distributions to the amount of 
levels included and to the possibility of having a prolate 
or oblate shape. 

Keeping just this quadrupole deformation we have 
compared the results of the coupled-channel calculation 
with the results of the frozen approximation. We show 
how the frozen approximation slightly overestimates the 
total cross section below the Coulomb barrier. Neverthe¬ 
less, the barrier distribution from the coupled channel 
calculation and from the frozen approximation show an 
overall good agreement for the present nuclei. 

Finally, we analyze the effect of adding an octupole 


deformation or vibration. Even though the differences 
in the form factors are not large and the same levels are 
included, the fusion cross section is sensitive enough to 
change from one coupling scheme to the other. Both 
cases considered, 16 0+ 144 Ba and 16 0 + 224 Ra, show slight 
differences in the barrier distributions. This difference is 
larger in the 224 Ra case. However, the final possibility of 
disentangle both distributions will depend on the avail¬ 
able accuracy in a hypothetical future experiment and 
the strength of the octupole deformation. 

Even if this difference may be not large enough for a 
solid conclusion towards clarifying a static octupole de¬ 
formation, it is worth to say that, to our knowledge, this 
is the first time that octupole deformation and vibration 
fully coupled to a quadrupole deformation are included in 
a coupled-channel calculation of a fusion reaction. This 
should be of great help in general for the analysis of nu¬ 
clear fusion with nuclei from regions with a large octupole 
deformation. 

Since the search for a large EDM needs large static 
octupole deformations, it might be possible to further 
ensure its presence by analyzing the subbarrier fusion 
cross section of some of the candidates. The possibility of 
finding differences in the barrier distribution could make 
subbarrier fusion a valuable tool to support the presence 
of static octupole deformations. 

To sum up, we would like to stress that fusion reactions 
have been of great help on understanding the structure 
of heavy ions for many years. In this particular case, it 
can also help the community to find the perfect candidate 
for large EDM’s and, consequently, to test models beyond 
the standard model. 
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